Scaling of the conductance distribution near the Anderson transition 
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The scaUng hypothesis is the foundation of our understanding of the Anderson transition. It 
has long been realised that the conductance of a disordered system is a fluctuating quantity and 
often assumed, but never demonstrated in the literature, that the conductance distribution obeys a 
single parameter scaling law. We present a clear cut numerical demonstration of the scaling of the 
conductance distribution in the critical regime. 



I. INTRODUCTION 

The single parameter scaling hypothesis of Abrahams 
et at is the basis of our understanding of the Ander- 
son metal-insulator transition in disordered systemsJj In 
Ref. I it was proposed that the zero temperature conduc- 
tance G = (e^ /h)g, as measured by the "Thouless num- 
ber" , obeys a single parameter scaling law. However, the 
large sample to sample fluctuations in the conductance 
of disordered systems were not explicitly considered. In 
the critical and localised regimes the fluctuations are of 
the same order as the mean conductance. (The rela- 
tion of the mean conductance to the Thouless number 
is discussed in Ref. |2[) This led to suggestions that 
the scaling hypothesis should be reformulated in terms 
of the typical cejaductancep or perhaps the distribution 
of conductance.QB 

For a disordered system of size L in d — 2 -\- e dimen- 
sions, Altshuler et at estimated the cumulants c„ of the. 
conductance distribution using a field theoretic method.tl 
At the mobility edge they found that 



(i) = 



{L/ey 



n > no 



(1) 



If single parameter scaling holds, the only relevant length 
should be the correlation length ^, and the appearance 
of the mean free path £ in the expression for the higher 
cumulants is unexpected. (In the insulating regime ^ is 
the localisation length, while in the metallic regime it is 
the correlation length.) Shapiro reconstructed the criti- 
cal conductance distribution from ([l]) and showed that, 
appearances to the contrary, these cumulants are consi&. 
tent with a single parameter scaling of the distribution.Q 
To apply results for c? = 2 -I- e to three dimensions we 
must make a questionable extrapolation to e = 1. While 
the Anderson transition occurs at weak disorder kp£ 1 
when e <C 1, it occurs at strong disorder kp£ « 1 in 
three dimensions. (Here kp is the Fermi wave number.) 
Comparison of the distribution obtained by Shapiro with 



numerical results shows that the behaviour of the conduc- 
tance distribution at large g, and also the non-universal 
behavioutjOf the higher cumulants in (|l]), is qualitatively 
incorrect .Ijtl 

To try to overcome the limitation to weak disorder 
Cohen et al. used a Migdal-Kadanoff type real space 
renormalization scheme.ll3 They found that the scaling 
of the conductance distribution is described by two pa- 
rameters; only in the limit of weak disorder is single pa- 
rameter scaling recovered. However, and as pointed out 
by Cohen et al., the Migdal-Kadanoff scheme involves an 
uncontrolled approximation and some of the results ob- 
tained with it are known not to be correct. For example, 
in the metallic regime the predicted conductance fluctua- 
tions are too large and in disagreement with the theoret- 
ically well established and experimentally verified phe- 
nomena of universal conductance fluctuations.EirEj They 
concluded that while the Migdal-Kadanoff scheme may 
be exact for hierarchical lattices ^jp three dimensions it 
is primarily of pedagogical value .E3 

Recent work on Anderson localisation in one dimen- 
sion hasJaighlighted the importance of a second length 
scale IsX^ Deych et al. demonstrated the existence of a 
crossover between single parameter and two parameter 
scaling regimes dependent on the ratio of Is to the lo- 
calisation length ^. Single parameter scaling is observed 
when ^ > Zs, and two parameter scaling when ^ < Ig. 
The implications of this result for Anderson localisation 
in higher dimensions are not yet clear. 

Numerical studies of the Anderson model have demon- 
strated single parameter scaling of the localisalipn length 
of electrons in quasi-one dimensional systems,Ej and also 
of the mean resistance, mean conductance and typi- 
cal conduct|ajice near the Anderson transition in three 
dimensions .Ea Of themselves these studies do not rule out 
a two parameter scaling of the conductance distribution. 
For example, in the Migdal-Kadanoff scheme the mean 
of the logarithm of resistance obeys a single parameter 
scaling law while at the same time the distribution obeys 
a two parameter scaling law. However, numerical stud- 
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ies have also demonstrated the existence of a universal 
size independent critical distributioii, that is accessible 
by varying only a single parameter.13 This observation 
makes a two parameter scaling unlikely; two parameter 
scaling should require that both parameters be varied si- 
multaneously to access the critical point. Nevertheless, 
given that the scaling hypothesis is the foundation upon 
which our understanding of transport in disordered sys- 
tems rests, we feel that a clear cut demonstration of the 
scaling of the distribution of conductance, free of ques- 
tionable extrapolations or uncontrolled approximations, 
is called for. 



II. METHOD 

Following Ref. ^ a single parameter scaling law for the 
conductance distribution Phid) of a three dimensional 
system of linear dimension L can be mathematically for- 
mulated as follows 



(2) 



where X is a parameter which must obey a single param- 
eter scaling law 



dlnX 
A\tlL 



= (3{X). 



(3) 



A limiting process is implicit in (^; we refer the reader 
to Ref. II for a detailed discussion. The parameter X 
need not be one of the moments of the distribution. 

At first sight it appears that we must know the func- 
tional form of the function _F in in order to verify 
single parameter scaling of the distribution numerically. 
In fact, this is not so. The procedure which we have 
adopted is to analyse the scaling of the percentiles of the 
distribution. The precise definition of the percentile gq is 



PL{g)Ag 



(4) 



where < g < 1. By establishing single parameter scal- 
ing for a representative set of percentiles we indirectly 
establish (|^) and (||), provided that the scaling of dif- 
ferent percentiles are consistent. When considering the 
percentiles it is not necessary to distinguish \ng or 1/ g 
as it is when considering average quantities. 

We have analysed the conductance distribution of the 
Anderson model numerically. The motion of the electrons 
is described by 



(5) 



<ij> 



where Cj{Ci) is the creation (annihilation) operator of 
an electron at the site i of a three dimensional cubic 
lattice. The amplitude of the random potential at site i 
is Wi. Hopping is restricted to nearest neighbours and 



its amplitude was taken as the unit of energy, V = 1. 
We assumed a box distribution with each Wi uniformly 
distributed on the interval [—W/2, W/2]. In what follows 
we refer to the strength of the potential fluctuations W as 
the disorder. The numerical method used is described in 
Ref. The two terminal zero temperature conductance 
was evaluated using the Landauer formula 



g 



2tTth 



(6) 



where t is the transmission matrix describing the.propa- 
gation of electrons from one contact to the other .Eao 

The conductance distribution depends on the system 
size L, the disorder W, the Fermi energy Ep and the 
boundary conditions. We set Ep — 0.5 and imposed 
fixed boundary conditions in the transverse directions. 
We accumulated data for the disorder range 15 < < 
18 and system sizes 6 < L < 18. At the extremes of 
the disorder range the localisation (coixelation) length is 
of the same order as the system size,B2l so that our data 
covers the critical regime. 

To estimate gq we simulated 1,000,000 realisations of 
the random potential and calculated the conductance for 
each realisation. (For L = 18 the number of realisa- 
tions was approximately 500,000.) We sorted the data 
into ascending order and our estimate of gq is then the 
n = [qA'^rfJth datum in this list where [x] is the integer 
part of X. When fitting the numerical data it is also 
necessary to have an estimate of the accuracy of the per- 
centiles. Following the standard method we used the bi- 
nomial distribution to estimate the likely accuracy of the 
percentile. We define 



An= y'Ndqil-q), 



(7) 



and look up the (n -I- An)th and (n — An + l)th data 
in the list and calculate the differences with gq. Our 
estimate of the accuracy is then the largest of these two 
differences. In practice, we found that the accuracy of 
all the percentiles were comparable, being of the order 
of 0.2%. The data were then fitted with the finite size 
scaling forms below by minimizing the statistic in the 
usual way. 

To fit the system size and disorder dependence of the 
percentile we supposed a single parameter scaling law 
but allowed for deviations from scaling due to an irrel- 
evant scaling variable and non-linearity of the scaling 
variables .113 We fitted the data to 



Ingq = F{tp,(j)) , 



(8) 



where ip is the relevant scaling variable and cf) is the ir- 
relevant scaling variable. We approximated this scaling 
function by its first order expansion in the irrelevant scal- 
ing variable 

\iigq^Foiij) + <f>F,i^). (9) 
We expanded each scaling function as a power series 
Fo (x) = In {gq)^ + x + a2X^ + ■■■ + a„„a;"'' (10) 
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Fi (x) = 1 + bix ■ 



(11) 



Here {gq)^ is the critical value of the percentile. The scal- 
ing variables were approximated by expansions in terms 
of the dimensionless disorder 



w={Wc- W)/Wc 



(12) 



where Wc is the critical disorder separating the insulating 
and metallic phases, 



(13) 



and 



(j)^Ly ((/-o + C^lW + ^2W^ + ■■■ + (l>n,w"*) . (14) 

The critical exponent v describes the divergence of the 
localisation (correlation) length as the transition is ap- 
proached. 



^ = ^± {fplW + 1p2W H ^i^n^W 



(15) 



The constants and hence the absolute scale of the 
localisation (correlation) length ^, cannot be determined 
from the fit. The decay of the irrelevant scaling variable 
with system size is described by the exponent y < 0. 
Redundancy in the definition of the fitting parameters 
between the coefficients in the expansions of Fq and Fi 
and the expansions of ijj and (p are eliminated by setting 
some of the expansion coefficients of Fq and Fi to unity 
as shown. This choice is also necessary if Fq and Fi 
are to be universal. The total number of parameters is 
Np = no + ni + + + 4. 

The minimum of was found using the DRNLIN rou- 
tine of the IMSL numerical library. The starting values 
of the fitting parameters supplied to DRNLIN are in the 
region v « 1.6, Wc w 16.5, y « —3, -01 « 1. We set 
the initial value of In (gq)^ to a value close to its best fit 
value by visual inspection of the raw data and all other 
parameters were initially zero. The results of the fitting 
procedure are not especially sensitive to the choice of the 
starting values. A number of fits corresponding to dif- 
ferent choices of tiq, ni, and are possible and a 
selection criterion is necessary. We set a cut off for the 
goodness of fit probability Q at Q = 0.1 and searched 
for the fit which requires the fewest parameters which 
satisfies this. Broadly speaking all sensible choices of noi 
ni, ny, and lead to consistent estimates of the criti- 
cal parameters. The goodness of fit and the accuracy of 
the fitted parameters were estimated using Monte Carlo 
simulations of synthetic data sets.Lil 



III. RESULTS 

Results for the q ^ 0.025, 0.17, 0.5 and 0.83 percentiles 
are shown in Table |. (For a normal distribution with 
mean fi and variance cr^ , these choices would correspond 
to the points — 2tT, n — a, fi, ii + a in the distribution.) 



Precise details of the fits are given in Table ||. The esti- 
mates of the irrelevant exponent are consistent with those 
obtained in Ref. |6| and are not shown again here. 

Data for the g = 0.17 percentile are plotted in Figure 
1. In Figure 2 the same data, after subtraction of cor- 
rections to scaling, are re-plotted to demonstrate single 
parameter scaling. This is done by plotting the corrected 
data as a function of the ratio of the systems size L to 
the localisation (correlation) length When displayed in 
this way the data fall on two different curves correspond- 
ing to the localised (lower curve) and the delocaliscd (up- 
per curve) regimes. The two curves are described by two 
scaling functions F^ and F^ which are derived from _Fo, 



F+(a;) = ln(.g,)^ + xi/'' + ...+ 



(16) 



F_ (x) = In Cgg)^ - x^l- + ■ • • + (-1)"° a„„x"«/^ (17) 

Data on the metallic branch follow 

L 



In 5, = F+ 



^4 



while data on the the insulating branch follow 

L 



^^gq = F- 



(18) 



(19) 



For completeness some representative data for the me- 
dian (q = 0.5) and q — 0.87 percentiles appear in Figures 
3 and 4 respectively. 




FIG. 1: The q = 0.17 percentile of conductance distribution 
for a disordered L x L x L systems versus system size L for 
disorder W in the range [15, 18]. The lines are the fit of Eq. 



TABLE I: The estimated critical values for each percentile 
and 95% confidence intervals. 



q 














0.025 


16.48 ± 


.02 


-4.14 ± .03 


1.56 ±.03 


0.17 


16.47 ± 


.01 


-2.55 ± .01 


1.56 ± .01 


0.5 


16.48 ± 


.04 


-1.08 ± .03 


1.59 ± .03 


0.83 


16.46 ± 


.05 


0.13 ± .03 


1.60 ± .04 


TABLE IL The details of the fit. Na is the number of data. 


q 


no n^, 


ni 


Np 


Nd 




Q 


0.025 


3 2 


1 


10 


179 


159 


0.7 


0.17 


2 2 


1 


9 


179 


155 


0.8 


0.5 


3 2 


1 


10 


179 


172 


0.4 


0.83 


3 2 


1 


10 


179 


180 


0.3 



We also analysed the q = 0.975 percentile of the dis- 
tribution but were unable to convincingly fit its systems 
size and disorder dependence. The origin of the difficul- 
ties may be the large corrections to scaling encountered 
for this percentile. Larger systems sizes will probably be 
needed for a definitive analysis of the high conductance 
tail of the distribution. 




10 0.01 fi.y 1 10 



FIG. 2: The data of Figure 1, after subtraction of corrections 
to scaling, re-plotted as function of the ratio L/^ to display 
single parameter scaling. The lines are the scaling functions 
(Q and described in the text. 

For the percentiles analysed the estimates of the criti- 
cal disorder and the critical exponent obtained from the 
scaling of different percentiles are consistent as required. 
The estimates of the critical exponent in Table □ are con- 
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FIG. 3: Data for the median conductance (g = 0.5 percentile), 
together with the best fit of Eq. ^, as a function of disorder 
W for systems sizes in the range I/=6-18. 



sistent with our previous estimates based on-the scaling 
of the localisation length in ouasi-lc? systemsjij scaling of 
higher Lyapunov exponents, E3lJ and scaling of the meaa 
conductance, mean resistance and typical conductance.t^l 
The estimates are also consisterrt, with numerical esti- 
mates reported by other authors. a'E3E3 



IV. CONCLUSION 

Our numerical results demonstrate single parameter 
scaling of the zero temperature conductance distribution 
in the critical regime of the Anderson transition in three 
dimensions. This result complements a previous demon- 
stration of the scaling of the mean-conductance, typical 
conductance and mean resistance.li3 

A two parameter scaling of the conductance distribu- 
tion, similar to that found by by Deych et al. for one 
dimensional systems might be recovered in the strongly 
localised regime. The localisation length diverges at the 
critical point while Ig, which is related to the integrated 
density of states, is always finite. Thus on approaching 
the critical point we should always find ^ > Zs and sin- 
gle parameter scaling should be observed. Far from the 
critical point, if ^ becomes less than Is, a two parameter 
scaling might appear. It remains to be seen, however, if 
the results of Deych et al. carry over to higher dimen- 
sions. 
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FIG. 4: The data for the q = 0.83 percentile after corrections 
to scaling are subtracted and plotted as function of the ratio 
L/^ to make single parameter scaling evident. The lines are 
the scaling functions (^|)and (^) described in the text. 



We would like to thank the Institute for Solid State 
Physics of the University of Tokyo for the use of their 
computer facilities. PM would like to thank the Japan 
Society for the Promotion of Science and Sophia Univer- 
sity for their hospitality and financial support, and is also 
grateful for support under APVT Grant No. 51-021602. 



E. Abrahams, P. W. Anderson, D. C. Licciardello, and 
T. V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 
D. Braun, E. Hofstetter, A. MacKinnon, and G. Montam- 
baux, Phys. Rev. B 55, 7557 (1997). 

P. W. Anderson, D. J. Thouless, E. Abrahams, and D. S. 

Fisher, Phys. Rev. Lett. 22, 3519 (1980). 

B. Shapiro, Phys. Rev. B 34, 4394 (1986). 

B. Shapiro, Phil. Mag. B 56, 1031 (1987). 

B. L. Altshuler, V. E. Kravtsov, and I. V. Lerner, Phys. 

Lett. A 134, 488 (1989). 

B. Shapiro, Phys. Rev. Lett. 65, 1510 (1990). 

K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 78, 4083 (1997). 

P. Markos, Phys. Rev. Lett. 83, 588 (1999). 

A. Cohen, Y. Roth, and B. Shapiro, Phys. Rev. B 38, 
12125 (1988). 

B. L. Altshuler, JETP Lett. 41, 648 (1985). 

P. A. Lee and A. D. Stone, Phys. Rev. Lett. 55, 1622 
(1985). 

A. Cohen, Phys. Rev. B 51, 10406 (1995). 

L. I. Deych, A. A. Lisyansky, and B. L. Altshuler, Phys. 



Rev. Lett. 84, 2678 (2000). 

K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 82, 382 (1999). 
K. Slevin, P. Markos, and T. Ohtsuki, Phys. Rev. Lett. 86, 
3594 (2001). 

J. B. Pendry, A. MacKinnon, and P. J. Roberts, Proc. R. 
Soc. Loud. A 437, 67 (1992). 

E. N. Economou and C. M. Soukoulis, Phys. Rev. Lett. 46, 
618 (1981). 

D. S. Fisher and P. A. Lee, Phys. Rev. B 23, 6851 (1981). 
A. MacKinnon and B. Kramer, Z. Phys. B 53, 1 (1983). 
W. H. Press, A. A. Teukolsky, W. T. Vetterling, and B. P. 
Flannery, Numerical Recipes in Fortran (Cambridge Uni- 
versity Press, Cambridge, 1992). 

K. Slevin and T. Ohtsuki, Phys. Rev. B 63, 45108 (2001). 
P. Markos, J. Phys. A: Math, and Gen 33, L393 (2000). 
A. MacKinnon, J. Phys. Cond. Matt. 6, 2511 (1994). 

F. Milde, R. A. Romer, M. Schreiber, and V. Uski, Eur. 
Phys. J. B 15, 685 (2000). 

S. L. A. de Queiroz, Phys. Rev. B 63, 214202 (2001). 



